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ABSTRACT 

We present threshold enhanced QCD corrections to inclusive processes such as Deep inelastic scat- 
tering, Drell-Yan process and Higgs productions through gluon fusion and bottom quark annihila- 
tion processes using the resummed cross sections. The resummed cross sections are derived using 
renormalisation group invariance and mass factorisation theorem that these hard scattering cross 
sections satisfy and Sudakov resummation of QCD amplitudes. We show how these higher order 
threshold QCD corrections improve the theoretical predictions for the Higgs production through 
gluon fusion at hadron colliders. 



Perturbative Quantum Chromodynamics(pQCD) provides a framework to successfully com- 
pute various observables involving hadrons at high energies. Wealth of precise data from hadronic 
experiments such as deep inelastic scattering(DIS) experiments at HERA complimented with the 
theoretical advances that lead to enormous success in computing the relevant observables very pre- 
cisely has given us a better understanding of the structure of hadrons and also the strong interaction 
dynamics in terms of their constituents such as quarks and gluons at wide range of energies. With 
these successes we can now predict most of the important observables with less theoretical uncer- 
tainty for the physics studies at present collider Tevatron in Fermi-Lab as well as at the upcoming 
Large Hadron Collider(LHC) in CERN [1]. 

The Drell-Yan(DY) production of di-leptons is one of the important processes at hadronic 
colliders. It will serve not only as a luminosity monitor but also provide vital information on 
physics beyond Standard Model(SM). The other process which is equally important is Higgs bo- 
son production at such colliders because it will establish the Standard Model as well as beyond SM 
Higgs [2, 3]. In pQCD, the DY production of di-leptons and Higgs boson production are known 
upto next to next to leading order(NNLO) level in QCD [4-16]. Beyond NLO, the Higgs produc- 
tion cross sections are known only in the large top quark mass limit. Apart from these fixed order 
results, the resummation programs for the threshold corrections to both DY and Higgs productions 
have also been very successful [17, 18](see also [19]). For next to next to leading logarithmic 
(NNLL) resummation, see [20, 21]. Due to several important results at three loop level that are 
available in recent times [22-27], the resummation upto N^LL has also become reality [28-31]. 
These results in fixed order as well as resummed calculations unravel the interesting structures in 
the perturbative results (for example: [31-34]). 

In the QCD improved parton model, the infra-red safe observables, say, hadronic cross sec- 
tions can be expandable in terms of perturbatively calculable partonic cross sections appropriately 
convoluted with non-perturbative operator matrix elements known as parton distribution func- 
tions(PDF). This is possible due to the factorisation property that certain hard scattering cross 
sections satisfy. The partonic cross sections are calculable in powers of the strong coupling con- 
stant gs within the framework of perturbative QCD because the coupling constant becomes small 
at high energies. On the other hand we only know the perturbative scale evolution of the parton 
distribution functions which otherwise are known/extracted from experiments. 

The fixed order QCD predictions have limitations in applicability due to the presence of various 
logarithms of kinematical origins. These logarithms become large in some kinematical regions 
which otherwise can be probed by the experiments. In such regions, the applicability of fixed 
order perturbative results becomes questionable due to the missing higher order corrections that 
are hard to compute with the present day techniques. The alternate approach to probe these regions 
is to resum these logarithms in a closed form. Such an approach of resuming a class of large 
logarithms supplemented with fixed order results can almost cover the entire kinematic region 
of the phase space. In addition, these threshold corrections are further enhanced when the the 
flux of the incoming partons become large in those regions. In the case of Higgs production 
through gluon fusion, the gluon flux at small partonic energies becomes large improving the role 
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of threshold corrections. Here,we consider the inclusive cross sections of hadronic cross reactions 
such as deep inelastic scattering, DY, Higgs production through gluon fusion and bottom quark 
annihilation and study the effects of soft gluons that originate in the threshold region of the phase 
space. In these processes, large logarithms are generated when the gluons that are emitted from 
the incoming/outgoing partons become soft. 

In [31], we extracted the soft distribution functions of Drell-Yan and Higgs production cross 
sections in perturbative QCD and showed that they are maximally non-abelien. That is, we found 
that the soft distribution function of Higgs production can be got entirely from the DY process 
by a simple multiplication of the colour factor Ca/Cf. In this article we extend the applicability 
of this method to the other important process, namely Higgs production through bottom quark 
annihilation. Using the soft distribution functions extracted from DY, and the form factor of the 
Yukawa coupling of Higgs to bottom quarks, we can now predict soft plus virtual part of the Higgs 
production through bottom quark annihilation beyond NNLO with the same accuracy that DY pro- 
cess and the gluon fusion to Higgs are known [28, 31]. This generalises our earlier approach to 
include any infrared safe inclusive cross section. In [31] we determined the threshold exponents 
D\ upto three loop level for DY and Higgs productions using our resummed soft distribution func- 
tions. In this present work we provide all order proof which establishes the relation between soft 
distribution functions and the threshold resummation exponents and demonstrate the usefulness 
of this approach to derive higher order threshold enhanced corrections for any infrared safe in- 
clusive cross sections. The results on DIS in this article provide a consistency check on various 
approaches (see [28, 31]) in the literature to study the soft part of the cross sections because some 
of the constants that go into our study were extracted from the deep inelastic scattering coefficient 
functions [25]. In the following, we systematically formulate a framework to resum the dominant 
soft gluon contributions in z space to the inclusive cross sections. Here the variable z is the ap- 
propriate scaling variable that enters the cross sections. To achieve this, we use renormalisation 
group(RG) invariance, mass factorisation and Sudakov resummation of QCD amplitudes as the 
guiding principles. Using the resummed results in z space, we predict soft plus virtual part of the 
dominant partonic cross sections beyond N^LO. The soft plus virtual corrections are also called 
threshold corrections. At the end, we show how these results affect the total cross section and 
improve the scale uncertainty for the Higgs production through gluon fusion. 

Since we are only interested in the effect of soft gluons, the infra-red safe observable can 
be obtained by adding soft part of the cross sections with the virtual contributions and perform- 
ing mass factorisation using appropriate counter terms. We call this infra-red safe combination 
a "soft plus virtual"(5v) part of the cross section. The soft plus virtual part of the cross section 
{A.^] p{z,q^,l^,l^)) after mass factorisation is found to be 



where ^^(z, e) is a finite distribution. The subscript P = 5 for Drell-Yan(DY) and Higgs 

productions and P = SJ for deep inelastic scattering. The symbol S stands for "soft" and SJ stands 
for "soft plus jet". For DY and DIS, / = ^(quark/anti-quark) and for Higgs production through 
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gluon fusion, / = g(gluon) and for bottom quark annihilation to Higgs boson, / = Z7(bottom quark). 
Here ^p(z, q^,i^,n^,z) is computed in 4 + e dimensions. 

%{z,q^,liWF,z) = (^ln(z^(4,//i//^e))' + ln|F^(4,!2',//',e)|')5(l-z) 

+2^|,{ds,q^,l?,z,^)-lmC\nTn{ds,n^,nl,z,^), I = q,g,b(2) 

In the above m = 1 for DY and Higgs productions and m = 1/2 for DIS. The symbol "C" means 
convolution. For example, C acting on an exponential of a function f(z) has the following expan- 
sion: 

Cefi^) = 5(1 -z) + ^f{z) + ^fiz) + ^/(z) + • • • (3) 

In the rest of the paper, the function /(z) is a distribution of the kind 5(1 —z) and 2),, where 

"In'(l-z)" 
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/ = 0,1,--- (4) 



and the symbol (8) means the Mellin convolution. Since we are only interested in the soft plus 
virtual part of the cross sections, we drop all the regular functions that result from various convo- 
lutions. F\ds, Q^jliP') are the form factors that enter in the Drell-Yan(/ = q) production, Higgs(/ = 
g^b) production and DIS(/ = q/q) cross sections. In the form factors, = —q^. For the DY, 
= is the invariant mass of the final state di-lepton pair and for the Higgs production, 

q^ = mjj, where mn is the mass of the Higgs boson. For DIS, = q^/z* ^^e virtuality of the 
photon or Z. The variable z in DIS is Bjorken scaling variable. In DY and Higgs production, 
z is the ratio of q^ over s, where s is the center of mass of the partonic system. The functions 
(^p{ds,q^,ILp',z) are called the soft distribution functions. The unrenormaUsed(bare) strong cou- 
pling constant dg is defined as 
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where gs is the strong coupling constant which is dimensionless in n = 4 + £, with n being the 
number of space time dimensions. The scale // comes from the dimensional regularisation in order 
to make the bare coupling constant gs dimensionless in n dimensions. 

The bare coupling constant ds is related to renormalised one by the following relation: 

Szds = Z(jjl)as{/jl) (6) 

The scale jjr is the renormalisation scale at which the renormalised strong coupling constant asijj^ 
is defined. The factorisation scale hf is due to mass factorisation. 

5e = exp{|[Y£-ln47i]} (7) 
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is the spherical factor characteristic of n-dimensional regularisation. 
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Z(A/^) = l+a,(A/^)^ + a^(A/^ 
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The renormalisation constant Z(^^) relates the bare coupling constant ds to the renormalised one 
asiiJ^) through the eqn.©. 

The coefficients (3o,|3i and (32 are 
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where the color factors for SU (N) QCD are given by 
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and nf is the number of active flavours. In the case of the Higgs production, the number of active 
flavours is five because the top degree of freedom is integrated out in the large rutop limit. 

The factors Z' {ds,^i^,fj^,£) are the overall operator renormalisation constants. For the vector 
current Z'i{ds,i^,i^) = 1, but the gluon operator gets overall renormalisation [35] given by 
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and for the bottom quark coupled to Higgs, we have 
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where the anomalous dimensions 7^ can be obtained from the quark mass anomalous dimensions 
[36] 
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The bare form factors F'{ds,Q^,f/,e) (without overall renormalisation) of both fermionic and 
gluonic operators satisfy the following differential equation that follows from the gauge as well as 
renormalisation group invariances [37^0]. In dimensional regularisation, 

.2 \ / n2 .2 
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where contains all the poles in 8. On the other hand, collects rest of the terms that are finite 
as 8 becomes zero. Renormalisation group invariance (RG) of the F'{as,Q^,iJ^,E) leads 



"^Mr \ Mr M J 



(15) 



The quantities are the standard cusp anomalous dimensions and they are expanded in powers of 
renormalised strong coupling constant ^^(/il) as 



(16) 



The total derivative is given by 



=^^^^ + -^^^77-21 (17) 



The RG equation for G' can also be solved and the solution is found to be 
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The finite function G\as{^\ l,e) can also be expanded in powers of asiQ^) as 

oo 

G^(a,(e2), i,e) = Y^^liQ^) ^/(e) (19) 

The solution to the eqn. (fT?l) can be obtained as a series expansion in the bare coupling constant in 
dimensional regularisation. The formal solution upto four loop level can be found in [24, 31] 

The coefficients G/(£) can be found for both I = q and I = g in [25] to the required accuracy 
in £. We have extended this to the form factor corresponding to Yukawa interaction of Higgs boson 
to the bottom quarks. We have presented the logarithm of this form factor, InF^ in the appendix 
A. All these form factors satisfy 



= 2(B{-8j,,^o-8,,bi)+fi+t^'s;'' 

V / k=i 

G/(£) = 2f5i-2 5/.,pi-5/,,y[)+/i-2po5/-^ + f;eV'' 

V / k=i 

G/(e) = 2 (^Bi - 3 5;,, ^2 - 6 i,t ii)+fi- 2|3ig/'^ - 2|3o (g/'^ + 2|3og/'') + £ e^'' 
G/(£) = 2(^5i-45/,,|33-6/,,7^^+/i-2|32^/'i-2(3i(^/'i+4pog/'2) 

-2po(^/'' +2|3og/'V4|32g/'3j + leV'' (20) 



Notice that the single poles of the form factors contain the combination 



2(^5/-6/,,/|3,_i-5/,,7f_i)+// 



at every order in a.,. The terms proportional — 2(5/,g / P,_i +5/ ;,7* J come from the large mo- 
mentum region of the loop integrals that are giving ultraviolet divergences. The poles containing 
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them will go away when the form factors undergo overall operator UV renormalisation through 
the renormalisation constants Z ^ which satisfy RG equations 

I oo 

^\-^\nZ\dsA.^\€) = £aUA/^)7ti (21) 

where £ ^ is set. The constants i [3, i and are anomalous dimensions of the renormalised 
form factors and respectively. The poles of the form factors containing terms proportional 
to ^.^j multiplied by |3/ are due to coupling constant renormalisation. It is now straightforward 
to apply this procedure to form factor of many operators in QCD that are of interest. After the 
overall operator renormalisation through Z ' and coupling constant renormalisation through Z, the 
remaining poles will contain only B/ and fj' in addition to the standard cusp anomalous dimen- 
sions a/. The constants are known upto order thanks to the recent computation of three 
loop anomalous dimensions/splitting functions [22,23]. They are found to be flavour independent, 
that is 5? = 5f . The constants fj are analogous to the cusp anomalous dimensions A' that enter 
the form factors. The cusp anomalous dimensions are flavour independent, = A^. It was first 
noticed in [41] that the single pole (in £) of the logarithm of form factors upto two loop level (a^) 
can be predicted due the presence of constants // because these f^' are found to be maximally 
non-abelien obeying the relation 

Jf = ft = ^ff (22) 

similar to A'-. In [25], this relation has been found to hold even at the three loop level. From 
the experience upto three loop level, we can now predict all the poles of the form factor at every 
order in ds from these constants A',B',f^ , their anomalous dimensions and the finite parts of lower 
order(in a.,) contributions to the form factor. The last equality in eqn. lOUb has been expanded again 
in terms of and g''^ with i,k = 1,2,3 so that the single pole can again be predicted. One is 
not sure whether this structure will go through beyond three loop until an explicit calculation is 
available. 

The coUinear singularities that arise due to massless partons are removed in MS scheme using 
the mass factorisation kernel r{z,/Jp,£) as shown in the eqn.Q. We have suppressed the de- 
pendence on ds and /j^ in P. The kernel r(z,;U^,£) satisfies the following renormalisation group 
equation: 

^l-^r{z,nl,e) = l:P{z,^il) ®r(z,/i2,e) (23) 

The P{z,^^f) are well known Altarelli-Parisi splitting functions(matrix valued) known upto three 
loop level [22,23]: 

00 

PizA) = l^alii4)P^'-'Hz) (24) 

i=l 
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The diagonal terms of splitting functions P^'^ (z) have the following structure 



5f+i6(l-z)+Af+i©o 
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where P^'J^ jj are regular when the argument takes the kinematic limit(here z ^ 1). The RG equation 
of the kernel can be solved in dimensional regularisation in powers of strong coupling constant. 
Since we are interested only in the soft plus virtual part of the cross section, only the diagonal parts 
of the kernels contribute. In the MS scheme, the kernel contains only poles in £. The kernel can be 
expanded in powers of bare coupling ds as 



.£ 

nzA^e) = 6(1 -z) + £4 i^^j 5^r«(z,£) 



(26) 



The constants T^'^ (z, e) expanded in negative powers of e upto four loop level can be found in [31]. 
The Tu{as,iJi]^,iJi^, z, £) in the eqn.© is the diagonal element of r(z,/i^, £) . 

The fact that l:^^ p are finite in the limit e ^ implies that the soft distribution functions have 
pole structure in £ similar to that of and F//. Hence it is natural to expect that the soft distribu- 
tion functions also satisfy Sudakov type differential equation that the form factors satisfy(see 
eqn.dni)): 



Kp \ as,j^,z,£ j+(~rp ^ai.,--2-,-^,z,e 



(27) 



where again the constants Kp contain all the singular terms and g/s are finite functions of £. Also, 
(^p{ds,q^,fj^,z) satisfy the renormalisation group equation: 



$/(a,,^^A'^z,e) = 



(28) 



This renormalisation group invariance leads to 



dpi V A' 



-A'{as{^ilmi-z) 
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z.e\=A'{a,{^im\-z) 



(29) 



If ^p{ds,q ,fr,z,£) have to contain the right poles to cancel the poles coming from F ,Z and F// 
in order to make A^J p finite, then the constants A ' have to satisfy 



(30) 
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Using the above relation, the solution to RG equation for Gp ^a^, z, is found to be 



Gp (a,.(A/i),%,z,e 



g/ {a,{q^), l,z,£) -6(1 -z) |j ^A' {a,{lVR)) OD 



7, 

With these solutions, it is straightforward to solve the Sudakov differential equation: 



where 



1 



(33) 



The above solution depends on m so that or equivalently Afp is finite as e ^ 0(see eqn.©). 

The constants K are independent of P and are determined by expanding Kp in powers of 

bare coupling constant ds as 

£ 

,2 \ _°o /,,2\'2 

(34) 



£ 

2\ °°/2\'T 

4(a.,^,z,£j =6(l-z)£ai,(^^j 5^^^«(£) 



and solving RG equation for Kp ('a?, ^,z,£ ) . We obtain 



Z^'(^)(£) 



Z^'(^)(£) 
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The constants Gp^^ (e) are related to the finite function G^p{as{q^) , 1 , z, e) through the distributions 
5(1 -z) and 'Dj. Defining Q ^{z) through 

.£ 

£a;f^^ii^)'"5^G/'('\£) = ta\{q\l-zf-)^U^) (36) 

;=1 \ ^^ J (=1 

we find 

G/'^'^e) = ^(-2|3o?p'i(e))+?^(£) 

Gp^'\^) = ^(4p^?p'i(e))+^(-Pi?^(e)-4po?^(e))+?P3(e) 

^/'^'^(e) = ^(-8|33^^^i(£))+l(^poPi?p'i(£) + 12|32^py£)) 

+^(-^p2^A(e)-2pi?^(£)-6po?p3(e)^ +?p'4(e) (37) 

The z independent constants Q pi{z) are obtained by demanding the finiteness of A^j'p given in 
eqn.lHl). Without setting £ = in cqn.©, we expand A"'p as 

/^lp{z.q\^WF,z) = Y^a',ipl)A']f{z,q\ijl,pl,£) (38) 

i=0 

Now, using the above expansion and eqn.Q, we determine these constants by comparing the pole 
as well as non-pole terms of the form factors, mass factorisation kernels and coefficient functions 
A^fp expanded in powers of £ to the desired accuracy. The factor m appearing in the solution 
turns out to be 1(1/2) for the DY and Higgs productions(DIS). This choice is governed by the 
number of incoming partons that are responsible for the soft emissions. Since G/(£)s in the form 
factors are found to satisfy a specific structure in terms of B', f',^i and 7^ as given in eqn. lODb . we 
find that the constants Q pi{£) also satisfy similar looking expansion containing these constants. 



^A(e) = -(fi+B{8p,sj) + t^'^pf 
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(39) 



The above structure indicates that all the poles including the single pole of the soft distribution 
function can be predicted from that of the form factors, renormalisation constants and the mass 
factorisation kernels. This is possible because we have now better understanding [41] of the struc- 
ture of even the single pole terms of the form factors. Notice that the terms proportional in the 
above expansion are due to the coupling constant renormalisation. The coefficients of single poles 
are proportional to the cusp anomalous dimension and the combination —{f + B^bp^sj) where 
the constants /f and B\ are process independent. The e dependent terms in g p{€) can be obtained 
from the fixed order(in Us) computations of cross sections and the finite parts of the form factors. 
At the moment, we know g p-^{t) to all orders in £, g p2{^) to order E and g p^{'e) to order z^. The 

lowest order term g p-^ (e) is known to all orders in £ because it is straight forward to compute the 
fixed order soft contribution. On the other hand, it is technically hard to determine 8 dependent 
parts of soft cross sections beyond the lowest order a^. We find. 








(40) 



where C/ = Cp for I — q^b and Q — Ca for / = 
be maximally non-abelien. That is, they satisfy 
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Interestingly these constants g s ii^) turn out to 




(41) 



Similarly, for the DIS, the constants • are found to be 
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The threshold corrections dominate when the partonic scaling variable z approaches its kine- 
matic limit which is 1. They manifest in terms of the distributions 5(1 — z) and 2),. Since the 
hadronic cross sections are expressed in terms of convolutions of partonic cross sections and the 
parton distribution functions, it is more convenient to study the threshold enhanced corrections 
in Mellin space where z ^ 1 corresponds to large A^. In Mellin space, all the convolutions 
become ordinary products and the 5(1 —z) distribution becomes a constant and the distributions 
'Dj become functions of logarithm of A^. The threshold resummation in Mellin A^ space has been 
a successful approach thanks to several important works given in [17, 18,42,43]. We show in the 
following how the soft distribution function 4>p(fl^,^-2,/|2,z,e) captures all the features of the A'^ 
space resummation approach. The exponents of the z space resummed cross sections get contribu- 
tions from both form factor as well as the soft distribution functions. The form factor contributes to 
5( 1 — z) part and the soft distributions functions contribute to 5( 1 — z) as well as to the distributions 
-Dj. Using 




(43) 



we can express the soft distribution function (for any m) as 





\l-zlJ^^l 
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(44) 



where 




(45) 
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For m = 1, that is, for DY and Higgs production, we can easily identify G/ [gs {q^{l — z)^) ,£) 
with the threshold exponent (a^ [q^{ \ — z)^)): 

oo 



i=l 



= 2Gs{as{qH'^-zf),e) 
Comparing the above result with the eqn. lBHb . we find 



(46) 



£=0 



Df = 2(?//£ = 0) (47) 

For m = 1/2, that is, for DIS, we identify G^j [us (^^(1 — z)) ,£) with the threshold exponent 
B'Disi^siqH-^-z))): 

oo 

1=1 



= G^j{as{q\l-z)),e) 
Comparing the above result with the eqn.(l36b. we find 



(48) 

£=0 



B'Dis,i=gU^ = 0) (49) 

It is evident from the above derivation that the allowed value m = 1/2 for DIS does not allow 
an exponent similar to D'{as{q^{l — z)^)) that appears in DY and Higgs productions. It is also 
straightforward to show that the so called conspicuous relation found in [26] follows directly from 
eqns. (l39l47l49l) : 

\D^-Bls,=B^ + fqmSg) (50) 

^ j=o 

where the constants ((3) and 5^^ can be obtained from the eqns. (l39l40l42l) . 

Few important remarks about the eqn.(|44l) are in order. The third line in the eqn. (|44ll cancels 
against ©o part of the mass factorisation kernel. The second line contains the right poles in £ to 
cancel those coming from the form factor as well as the mass factorisation kernel. In addition, 
it contains the terms that are finite as £ becomes zero through the constants g p-{e). These finite 



constants contribute to soft part of the cross section at higher orders. Hence, adding the eqn. (|44ll 
with the renormalised form factor and the mass factorisation kernels, performing the coupling 
constant renormalisation and then finally taking Mellin moment, we reproduce the resummed result 
given in [17, 18,42,43] when £ ^ 0. 

In [31], it was observed that the soft distribution functions 4>/(a^,^^,;U^,z,£) of DY and Higgs 
production are maximally non-abelien and hence the entire soft contribution of Higgs production 
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can be predicted from that of DY and vice versa. In fact the £ dependent parts of Q g ex- 
tracted from DY are needed to predict the soft part of the Higgs production. In addition, using the 
resummed result given in eqn.lH]), and the available exponents ^/(e), ^ i (e), one can also predict 
part of the higher order soft plus virtual contributions to the cross sections. The available exponents 
are 
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in addition to the known P,; (z = 0, 1, 2,3), the constants in the splitting functions A,-, 5,- (z= 1,2,3), 
the maximally non-abelien constants (z = 1 , 2, 3) and the anomalous dimensions 7f(z' = 0, 1 , 2, 3). 
For I — q,g, the constants and gf-^ are known for j = 2,3 also (see [24]). 

Using the resummed expression given in eqn.lH]) and the known exponents, we present here 
the results for A^p'^ for DY, Higgs production and DIS. As we have already mentioned in the 
beginning, for DIS, it is just a consistency check because the three loop form factors were derived 
from the DIS coefficient functions. The Drell-Yan coefficients A^^'i'"* and Higgs coefficients A^\^'^ 

and A^^^^'^ are known upto NNLO (z = 0, 1,2) from the explicit computations (see [4]- [16]). For 

N^LO for I = q,g, a partial result A^i'^'^\ i.e., a result without 5(1 —z) part is available from 
the work of [28]. Using the universal behavior of the soft distribution function for P = 5, in the 
reference [31], the entire A^^^*^'^ for I = q,g upto NNLO as well as partial N^LO for Drell-Yan and 

Higgs through gluon fusion, that is the coefficients Af'j^^ for I = q,g were reproduced. This was 
achieved using the resummation formula given in eqn.lH]). In this article, we predict for the first 
time the N^LO contribution A^^^^'^^ to Higgs production through bottom quark annihilation with 

the same accuracy that the coefficients A^^^'^^'' for I = q,g are known. In addition, we extend this 
approach to N^LO order where we can predict partial soft plus virtual contribution coming from all 
'D j except J = 0, 1 for Drell-Yan N^LO coefficient A^^^^^\ gluon fusion to Higgs N^LO coefficient 

'^g 5^^ and bottom quark annihilation to Higgs boson N^LO coefficient A^^^^^-*. Like N^LO Aj^^*-^"*, 
here also we can not predict 6(1 — z) part. These results are presented in the Appendix B for 
lj^= juj^ = q^. Using our resummed formula and the DIS exponents g^j{e), we have reproduced^ 
the partial N^LO and N^LO coefficients A^^^^"*, A^^^j"* for DIS given in [26,45] and the lower order 
results (see [44] for NNLO). The computation of these soft plus virtual contributions to DY, Higgs 
productions and DIS involves convolutions of distribution functions "Dj with j = 0, ...,7. We have 
computed them using the following integral representation given in [46]. 

= lim lim fl-V^-V^ C dx C dy{\-xy^-\\-yf^-^ 
J e^o a,b^i \da J \db J e'+J Jo Jo 

'after fixing the typos in eqn.(4.19) of [45] and in eqns.(5.6,5.7) and eqn.(5.9) of [26] 
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X 1 6(z - xy) - 6(z -x)-d{z-y) + b{l-z)^ (51) 

The above integral reduces to a set of Hypergeometric functions and Euler Gamma functions which 
can be expanded around £ = to desired accuracy. We are interested only in those terms which 
are proportional to the distributions 5( 1 — z) and for our analysis. For this purpose, we have 
derived a formula to compute the convolutions of distributions for any arbitrary /, j using the above 
integral representation: 



, dVfdy 1 fa + b 
-Di^-Dj = hm km — — — — 



r{l+ae)r{\+be) _^ 



r{l + {a + b)e) 

k / \ k 

r 1 . ^ ({a + b)E] 

X 



1 °o ((a + b)e) 1 °o ((a + b)e' 
_i_5(l-z)+y lA-z), + 1^2), 



k=i 



^ (aef b ^ {bef , A 



a + b ,v. , . 



where R{a,b,z,£) is the remaining regular function. It is now straightforward to obtain A^^^^'^ for 
i= 1 , . . . , 4 for both Higgs(/ = g,b) and DY(7 = q) productions. 

The impact of partial soft plus virtual parts of N^LO and N^LO contributions to Higgs produc- 
tion through gluon fusion at LHC is presented in figure (1). The Higgs production cross section is 
given by 

U I C dy^ab{yA)^abi-.mWF.lA) (53) 
where x = mjj/S, N = 3 and the factor Gb can be found from [35]. The flux <i>afe(j,//^) is given by 

^abiyA) = — faiwA) fb (54) 

where fa{w,i/p) is the parton distribution function. The partonic cross section A^^ contains both 
soft plus virtual as well as hard contributions: 

^^biz,mWRA)=^7A^,mj,AA)+Kb'"''i^^^^^^ (55) 

We choose the center of mass energy to be \/S = 14 TeV for LHC. The standard model pa- 
rameters that enter our computation are the Fermi constant Gf = 4541.68, Z boson mass Mz — 
91.1876 GeV, top quark mass = 173.4 GeV. The strong coupling constant CXj(jU^) is evolved 
using 4-loop RG equations depending on the order in which the cross section is evaluated. We 
choose a^^(Mz) = 0.130, af^^(Mz) = 0.119, af^^^(Mz) = 0.115 and af^{Mz) = 0.114 for 
/ > 2. We use MRST 2001 LO for leading order, MRST2001 NLO for next to leading order and 
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LHC(14TeV) 

N**LO(pSV) 

N^LO(pSV) 
NNLO 
NLO 
LO 




1.4 



1.2 



0.8 




Figure 1 : Total cross section for Higgs production through gluon fusion at LHC, and its scale 
dependence for the Higgs boson of mass mn = 150 GeV with // = mr^iuq = hh- The abbreviation 
"pSV" means partial soft plus virtual. 



MRST 2002 NNLO for N'LO with / > 1 [47,48]. In the first plot of fig.(l), we have shown the 
cross sections in pb for various values of Higgs mass. For LO,NLO and NNLO we used the ex- 
act results which contain both soft plus virtual as well as regular hard contributions. For N'^LO 
(i=3,4), we use only soft plus virtual results extracted from the resummed formula. Here we have 
set Hp = hr = mn. We find that the inclusion of N'LO(i = 3,4) does not change the cross section 
much confirming the reliability of the perturbative approach. In the second plot of fig.(l), we have 
plotted the scale variation of the cross section using the ratio: 



m 



I) 



(56) 



It is clear from the second plot of fig.(l) that the inclusion of N'LO(i = 3,4) soft plus virtual 
contributions reduces the scale ambiguity further. 

To summarise, we have systematically studied the soft plus virtual correction to inclusive pro- 
cesses such as DIS, DY, Higgs productions through gluon fusion and bottom quark annihilation. 
Using renormalisation group invariance and Sudakov resummation of scattering amplitudes and 
the factorisation property of these hard scattering cross sections, the resummation of these cor- 
rections has been achieved. We have also shown how these resummed distributions are related to 
resummation exponents that appear in Mellin space. Using our resummed results we predict par- 
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tial soft plus virtual cross sections at N^LO and N^LO. We have also shown the phenomenological 
consequences of such results for Higgs production through gluon fusion at LHC. 

Acknowledgments: The author would like to thank Prakash Mathews for discussion. 



A Appendix 

In this appendix, we present the logarithm of the form factor InjF for the Yukawa interaction 
of Higgs boson with bottom quarks upto three loop level expanded in powers of 8 to the desired 
accuracy for our computation. 
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B Appendix 



In this appendix we present soft plus virtual parts of Drell-Yan production of di-leptons and Higgs 
productions through bottom quark annihilation and gluon fusion. 
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